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The investigation of the quantum systems having so-called non-Hermitian Hamiltonians 
has been a great interest because of its theoretical contributions to quantum mechanics [1- 
15]. The study of such quantum systems has been received many applications, especially in 
quantum field theories [16-18], and nuclear theory [19]. The solutions of the non-relativistic 
and relativistic equations with non-Hermitian Hamiltonians having real or complex energy 
spectra have also been studied by many authors by using different methods [3, 6, 11, 20]. 
Non-Hermitian Hamiltonians satisfy the condition given by OHO^ 1 = OHO = H, where the 
operator O = PT is an operator combined with parity, and time-reversal transformations, 
respectively. So they are related with PT-symmetry [8]. If a potential V(x) has the condition 
written as V(—x) = V*(x) under the transformation of x — > —x, and % — > —i, then it is said 
that the potential is a PT-symmetric. 

The relativistic wave equations, especially Dirac equation, give numerous important re- 
sults, and explanations in the view of quantum mechanics. One of them is the spin-orbit 
coupling, and related topics, such as the Hall effect, spin torque, Zitterbewegung, etc [21- 
23]. Many authors have also been studied the solution of the Dirac equation with different 
potentials in different theoretical backgrounds [24, 25]. 

Recently, the generalization of the solutions of the relativistic and non-relativistic wave 
equations with PT-symmetric potentials in the case of the constant mass to the case of 
the position-dependent mass become very attractive research topic. Many authors have 
studied the effects of the position-dependent mass on the solutions of the above equations 
[26-36]. The new formalism based on the position-dependent mass is a useful ground to 
study the properties of some physical systems, such as quantum dots [37], semiconductor 
heterostructures [38], quantum liquids [39]. 

In this letter, we intend to study the effects of the position-dependent mass on the solu- 
tions of the one-dimensional Dirac equation by using the map between the mass, and po- 
tential functions without the spin-effects [43]. Our aim, by taking two examples of the mass 
distributions, is to show that the two-component approach [40, 41] to the one-dimensional 
Dirac equation can be used in the position-dependent mass formalism. We obtain two new 
non- PT-symmetric, and non-Hermitian potential functions produced by using the map be- 
tween the mass, and potential functions. 
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The relativistic Dirac equation in the absence of an external potential can be written as 

(i^-rn)* = 0, (1) 

where m is the rest mass of the particle, and 7 M (/i = 0, 1, 2, 3) are gamma matrices (h — c — 
1). The Dirac equation for a particle moving in an external potential V(x) in one-dimension 
can be written as [40] 

[a .p + j3m(x) — (E — V(x))] ip(x) = , (2) 

where E is the relativistic energy of the particle, p is the momentum operator, m(x) is 
the mass function and a, [3 are 2x2 matrices which are set to Pauli matrices a 3 and ai, 
respectively. The Dirac wave function, ip(x), is decomposed into an upper ip + {x), and lower 
component ip~(x), so that 

#r)=(^V (3) 
V ip (x) J 

We have the following set of two coupled differential equations by inserting Eq. (3) into 
Eq. (2) 

di>+(x) 



dx 
dip~(x) 



+ [m(x) + (E- V{x))}^j-{x) = , (4) 
+ [m(x) -(E- V(x))]^ + (x) = . (5) 



dx 

We define the wave functions f(x) and g(x) in two-component approach which makes 
possible to obtain the solutions easily without the spin effects [40, 41] as 



g(x) ) V 1 — i J \ i])~{x) 



(6) 



and insert into Eqs. (4) and (5), we have 

df(x) 



dx 
dg(x) 



- i[E - V(x)]f(x) + im(x)g(x) = , (7) 
+ i[E - V(x)]g(x) - im(x)f(x) = . (8) 



dx 

By using the transformation written in terms of mass function m(x) in Eq. (7) 



f(x) = jM^j<f ) (x), (9) 
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one gets 



. d(f)(x) 
dx 



im (x) /ri T T . . . 
+ (E - V(x)) 



4>{x) - y/m(x) g{x) = . (10) 



2m(x) 

where prime denotes the derivative to the spatial coordinate. To obtain a Schrodinger-like 
equation for (f>(x), one sets the potential term V(x) in Eq. (10) as 



which gives 

d 2 (f)(x) 



V eff (x)<i>(x) = -E 2 <P(x) . (12) 



dx 2 

This is the Schrodinger equation having the energy eigenvalue E 2 under the effect of 
the potential V e ff(x) = m 2 (x). By choosing a suitable mass function, one can construct 
a potential V(x) from Eq. (11). By solving the Schrodinger-like equation given by Eq. 
(12), one obtains the energy spectra and corresponding spinors of the Dirac equation. The 
Schrodinger and Schrodinger-like equations can be solved by using the calculation procedure, 
which suggests a hypergemetric type equation in the following form [42] 

dz 2 a(z) dz a z (z) 

where cr(z), a(z) are polynomials at most second degree and f(z) is a first degree polynomial. 
By using the transformation ip(z) = x(z)y(z), one gets 

.«^ + r W ^+A !/W = 0, (14) 

The first part of the total wave function x{ z ) is defined as [42] 

X(z) cr(z) ' 

and the other part of the solution y(z) is given by the Rodrigues relation 



(15) 



a d n 

y n (z) = -^—[* n (z) P (z)}, (16) 

where a n is a normalization constant, and p(z) is the weight function obtained from the 
following relation [42] 

j- z [a{z)p{z)]=r{z)p{z). (17) 



The polynomial tt(z) in Eq. (15), and the parameter A required for the method are given 

by 



a (z) - r(z) 
7r{z) = ± 



'°JM^M) 2 - Hz) + k „ {2) . (18) 



\ = k + 7r'(z). (19) 

The discriminant of the expression under the square root in the polynomial 7r(z) in Eq. (18) 
must be zero, which defines the constant k. Thus, a new eigenvalue equation becomes 

A = A n = -nr'(z) - a"(z) . (20) 

where t(z) = f(z) + 2n(z) has a negative derivative. The energy eigenvalues are obtained 
from the Eqs. (19) and (20). 

We prefer the following mass function as a first case 

m{x) =m (l + e- 5x ), (21) 



which gives an exponential type, non-PT-symmetric, non-Hermitian potential 

iS e~ Sx 
2~I+e 



^) = ~r^- ( 22 ) 



By using the mass function, Eq. (12) gives 

fpp. + {e 2 -ml- 2m 2 e- Sx - m^- 2 ^}^^) = , (23) 
By inserting the new variable z~ l = e Sx , one obtains the following equation 

d 2 (p{z) 1 d<f>(z) 1 r~ 2 2 2 2 2\,( \ n / 0/1 n 

where P 2 = E 2 /5 2 , and a 2 = ml/5 2 . Comparing Eqs. (24) and Eq. (13), and using the 
parameters —a\ = —a 2 , — a 2 = —2a 2 , and — e = E 2 — a 2 , we obtain 

f(z) = 1 ; <j(z) = z ; a(z) = —a\z 2 — a 2 z — e . (25) 

To obtain the polynomial tt(z), one inserts these polynomials into Eq. (18), with o'(z) = 
1, then we get 

ir(z) = ±^/a lZ 2 + (k + a 2 )z + e , (26) 
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We obtain two solutions for k as k\ )2 = —a 2 ±2^aie , from the condition that the discriminant 
of the expression under the square root has to be zero. Eq. (18) gives 



n(z) 



±[ y/a^z + yje ], k — > ki , 
±[y/aTz — ^e], k^k 2 . 



We find the polynomial t(z) by choosing k 2 as 

t(z) = r (z) + 2n(z) = 1 + 2^ -2^z, (27) 

which gives r'(z) = — 2 v /a~T < 0. Thus, an energy eigenvalue equation is obtained from Eqs. 
(19) and (20) as 

{2n+ l) v V + 2 v /^F + a 2 = 0, (28) 



S 2 



2mn 



E = ±^m 2 - — [2n + l+ -^j . (29) 

It can be seen that the energy levels of particle and antiparticle are symmetric about 
zero. We observe that it is obtained zero energy in the case of the limit 5^0 which is the 
energy level of the following equation derived from Eq. (23) 

+ ^(x) = , (30) 

where u 2 = E 2 — itLq . The Dirac Hamiltonian with a general scalar potential, V'(x), gives 
always zero-energy solutions in the ultrarelativistic limit [43], where the upper and lower 
components given in Eq. (2) are written as 

ip-(x) ~ e +(mx+h(x)) _ ^ 

with h(x) = j x V'(y)dy , m denotes the rest mass of particle. The normalization conditions 
for zero-energy eigenstates are discussed in Ref. [43]. The solutions of Eq. (30) have the 
same form with the ones given in Eqs. (31) and (32) for the limit 5^0, i.e., h(x) = for 
m — > 2m . One has to find first the weight function p{z) to obtain the Dirac wave function 
4>{z). Using Eq. (17), we get 

p( z ) = e ~ 2 ^ z z 2 ^ 7 (33) 



and from Eq. (16), we obtain 

y n ( z ) = On r Wa 

V e -2^aTz Z 2^F dz n [ J V ; 

From the last equation, we can write y n (z) in terms of the generalized Laguerre polynomials 
as y n (z) ~ L^ 2 ^\z) [43]. The other part of the solution is obtained from Eq. (15) as 

X(z) = e-^ z z^ 7 . (35) 

Finally, we obtain the unnormalized wave function as 

(f) n (z) = e-^ z z^ L^\z) . (36) 

Thus, we write the upper component from Eq. (9) as 



f n (z) = ^Jm (l + z) e'^ z z^L^\z) (37) 
and the lower component from Eq. (8) as 

It is seen that the mass of the Dirac particle contributes to the eigenfunctions, and the 
upper component depends on the energy of the particle. 
We choose the second mass distribution 

™(*)=iq^> ( 39 ) 

which gives following exponential type, non-PT-symmetric and non-Hermitian potential 
from Eq. (11) as 

"<*> = *tt^- < 40 > 

It has the form of the Woods-Saxon potential for iS/2 — > — V [45]. Substituting Eq. (39) 
into Eq. (12), and using the new variable z — 1/(1 + e~ 5x ), we obtain 



where — P 2 = E 2 /5 2 . Comparing the Eqs. (41) and Eq. (13), we get 

f(z) = 1 - 2z ; <r(z) = z(l - 2) ; cr(*) = -£ 2 - a 2 z 2 . (42) 



Following the same procedure as in last section, we get the energy spectra 
1 



8m 2 - [V W 4m 2 ~ S(2n + l)] 2 - ( 4m ° )' . (43) 

Aj V ^5 2 + Ami - S(2n + 1) 1 



E = ±- 

4< 

and the corresponding wave functions are written in terms of the Jacobi polynomials 

Mz) = ^ (1 - z)' M Pi 2 *'- 2M \l - 2z) . (44) 

The lower and upper components g n (z) and f n {z) are written by using Eqs. (44) and (8) 
respectively as 

/„(*) = v^^ +1/2 (1 " zY M Pi 2 *>~ 2M Hl - 2z) , (45) 

and 

" 1 



9n(z) = z^ 2 (l-z)- M Pi 2 ^- 2M \l-2z){ iS 



m 



2 ;i-z) + e(i-z) 



+ Mz - z(l -*)(„ + 1 + 2E - 2M) gs.-.M, ' 



+ ^}. (46) 



where M = \/ E 2 + a 2 . 

The energy levels of the particles and antiparticles are symmetric about zero in Eq. (43). 
The non-PT-symmetric and non-Hermitian potential given by Eq. (40) has a real energy 
spectrum under the condition that 8ml > W $ 2 + ^ m o ~ &{2n + l)] 2 + 



v /<5 2 +4mg -<$(2n+l) 

The energy expression in Eq. (43) gives zero energy in the limit of 5 — > as in the first mass 
distribution case. Both of the eigenfunctions are dependent on the rest mass of the particle. 

In summary, we have investigated the position-dependent mass energy spectra, and the 
corresponding wave functions of the Dirac equation in two-component approach. We have 
transformed it into a Schrodinger-like equation with the help of a transformation applied on 
one of the Dirac wave functions by using two different mass distributions. Thus two new 
non-PT-symmetric and non-Hermitian complex potentials are produced. We have found 
that the Dirac equation can be turned into a Schrodinger-like equation by using a suitable 
transformation depending on the mass function m(x) in two-component approach. This 
formalism gives us two new non-PT-symmetric and non-Hermitian exponential and inverse- 
exponential potentials. We have obtained a complex energy spectrum for the exponentially 
potential given by Eq. (22), and the potential given by Eq. (43) has a real energy levels if 
iml > ly^Ti^ - S(2n + 1)P + ( Vsl+i ^ + J . 
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